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Abstract 

Filtration was produced in a small-scale physical model of a granular porous medium of cylindrical shape. The 
same volume flow was obtained either applying a difference in hydrostatic pressure or in osmotic pressure. In 
the first case a process of sustained erosion ending in a hydraulic short-circuit was observed, while in the second 
case the material remained stable. This paradoxical strength behaviour is explained using some results from 
differential geometry, classical field theory and thermo-kinetic theory. The fracture process of a continuous 
matrix in a porous medium under the combined effect of filtration and external mechanical loads is then 
considered. The obtained results can be applied to the flexural and compressive strength of wet concrete. 
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1. INTRODUCTION. 

Filtration was produced in a small-scale physical model of a granular porous medium 
of cylindrical shape, in the framework of an electro-osmotic experiment. 

A capillary glass tube (1 mm. of diameter) was filled with a matrix of densely packed 
grains (with linear dimensions between 0,1 um. y 1 um.). A cylindrical portion of 1 cm. of 
length was cutted from the tube, and without any aditional mechanical barriers at its ends, it 
was used as an horizontal filtration barrier between two water chambers (fig. 1). 

In order to determine the global hydraulic permeability of the filter, a difference in 
hydrostatic pressure was applied across it. The resulting volume flow was measured by the 
displacement of the air bubble in the long capillary tube also shown in fig. 1. 
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Fig. 1 The experimental setup 


Under the applied head of hydrostatic pressure, the material remained apparently 
stable during fifteen or twenty minutes. Then, the whole matrix was suddenly ejected towards 
the chamber that was receiving the volume flow: a true hydraulic short-circuit was produced 
between the two chambers. 

As it was not possible to measure accurately enough the volume flow by the bubble 
method to determine the desired hydraulic permeability, a second method had to be used. 

The same hydraulic permeability can be determined using an osmotic flow across the 
filter (Meschia and Setnikar, 1958). So, in order to produce a difference in osmotic pressure, 
the right chamber was filled with a 1 M sucrose solution, and the hydrostatic pressure at the 
tube level was mantained the same in both chambers. Now, for approximately the same 
volume flow, the granular material remained stable, not only during the whole experiment 
(three to four hours), but even after six or seven hours, until the system was disarmed in order 
to prepare it for a new experiment. It wasn’t possible to detect any sucrose in the left chamber, 
even after many hours of having stopped the osmotic flow by adding an excess of hydrostatic 
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pressure on the right chamber. As a consequence, the sucrose doesn’t enter in significant 
quantities in the bulk of the interconnected pore space of the filter, so that neither the strength 
of the bulk of the granular matrix nor the local flow characteristics can be modified by the 
sucrose. 

Then, the observed difference in strength behavior of the same matrix, depending of 
the force that produces the flow, seems to be somewhat paradoxical. This is so because. it is 
generally admitted that in the same porous material, with the same geometry and the same 
boundary conditions (in our case, free water at the ends of the filter and an impermeable glass 
boundary in its curved face), two equal global volume flows (one osmotic and the other 
hydraulic) must result from the same field of filtration velocities in the space of 
interconnected pores (Bear, 1972; de Marsily, 1986). 

But to the same field of filtration velocities it corresponds the same field of filtration 
pressures, that is, of forces per unit volume that act on the solid matrix during filtration, 
pushing in the flow direction. Then, under the same externally applied loads, the mechanical 
behavior of the bulk of the granular matrix should be the same, irrespective of the origin of 
the force that causes the flows. 


2. FILTRATION, EROSION AND OSMOSIS. 


However, a significant difference in strength behavior appears when sucrose is added 
to the right chamber. It seems reasonable to suppose that this difference is due to some 
process that has its origin in the boundary between the matrix and the water of the right 
chamber. This process can continue freely when the force that produces the flow is a 
difference in hydrostatic pressure, but it must be blocked when this force is a difference in 
osmotic pressure. 

Both conditions are satisfied assuming that an erosion tunnel grows, at first 
inadvertently, in the bulk of the matrix, from the water-matrix interface in the right chamber 
and against the flow, when the force is a difference in hydrostatic pressure between the 
chambers. This process is slow at its beginnings, so that the volume flow through the filter 
remains practically constant. But as the tunnel lengthens, its lengthening is produced faster 
and faster, and in the last part of the process the volume flow raises suddenly up to the 
observed hydraulic short-circuit. 

When the volume flow is due to a difference in osmotic pressure, the tunneling 
process isn’t produced because, as will be shown below, a mechanical force appears in the 
transition region between the sucrose solution and the bulk of the filter. This force counteracts 
the effects of the filtration pressure and blocks the beginning of the erosion process. 


3. MODELLING TOOLS. 


To model the aforesaid erosion process in the filter we will use three kinds of tools: 

(1) An average form of the equations that govern the transport processes in the pore space 
and the stresses in the solid matrix. 

(2) Some results for solenoidal and almost gradient fields, taken from classical field 
theory and differential geometry, and applied to the average field variables. 

(3) A thermodynamic argument and an approximate expression for the hydrostatic 
pressure of a liquid solution, as a function of the concentration of its components, that stems 
from the molecular theory of fluids. 


3.1 Average field equations. 
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In our case the flow can be considered as incompressible, steady and with negligible 
inertial effects in comparison with viscous forces (Stoke’s flow). 

In fact, the fields of velocity and stress in the fluid and the stress field in the solid 
matrix fluctuate randomly along distances of the order of the local pore’s diameter d (for a 
definition of the local pore diameter see Scheidegger, 1963). However it is possible to smooth 
out these fluctuations by averaging the field variables over regions of characteristic length/ , 
much greater than d and called Darcy’s scale. Selecting /, at its turn, much less than a 
characteristic length Lof the whole porous body, we obtain averages that can be considered 
as true local values (in our case we can estimate d = 10°’ m and L =~ 10° m, so that 
1 = JdL = 10° m will be suitable for averaging purposes). 

For each point in the porous medium we consider an averaging region that has this 
point as its center, and then we assign to the point the averaged values of all the field 
variables, with independence of the true location of the point (in porous space, in solid matrix 
or in the boundary between them). We thus obtain a vector field of filtration velocities V;, a 
scalar field of averaged pore hydrostatic pressure p, a scalar field of porosity ® and a tensor 
field of averaged matrix stresses GC. 

The emergent equations at Darcy scale are: 

Vev, =0 (v, solenoidal) [1] 


An 


Vp =- = eV Zz, (Darcy filtration law) [2] 
u 


where V e() represents the divergence, V() the gradient, ethe scalar product, 
T, =p + pgz is the filtration potential (f is water density, g is gravitational 
acceleration, z is the vertical position fron a reference horizontal plane, z increasing 
upwards), is water’s dynamic viscosity and K is a symmetric tensor with positive 


eigenvalues that express the local hydraulic conductivity of the porous medium. 
From [1] and [2] we obtain the following elliptic equation: 


Ve(KeVz,)=0 [3] 
At an impermeable boundary (as the glass-filter interface) V; « n = 0, ni being the 
unit normal vector to the interface at the considered point, so that 
fieKeVz, =0 [4] 
At a free water interface (as the left and right filter-water chamber interfaces) 7, must 
be constant: H, =C [5] 
The average force on matrix, per unit of volume of porous medium is: 
Pp; =-@®Vz, _ (filtration pressure) [6] 
Here @ is the (volume) porosity field, so that p, is everywhere orthogonal to a free water 
interface. The equation: 
Pp; =-Veo, [7] 
relates the filtration pressure with the average matrix stresses due to the filtration process, G;. 
The complete average stress tensor in the solid matrix verifies the emergent static 
equation 
Ve(G-®pl)=X [3] 
Here | is the unit tensor and X is the average external force per unit volume. 
Then, the force acting on an element of surface of area A A and unit normal vector n, 
at the scale of Darcy is given by the formula (Fliigge, 1972): 


AF=AA(6-pOleii [9] 
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3.2 Solenoidal and almost gradient fields. 


Let us consider a vector line of the fieldV,. If t, (s)is the unit tangent vector at the 
point of intrinsec coordinate s (s being the arc length), given any two points Q and P that 


belong to the line, we have the relation: 
sQ 
- Vet, (s) ds 


v:(Q)=v,@e * [10] 
(Herev, = el). Equation [10] is Bjorgum characterization for solenoidal fields (Ericksen, 
1963). If we assume that the filter is isotropic (although not necessarily homogeneous), then 


is “ K i ; 
K = KI reduces to a scalar field and Vv; = - — V a, is now an almost gradient field. 
u 


Now each tangent unit vector t, at a point of a filtration line is at the same time a unit 
normal vector to the surface of constant filtration potential that pass through the same point. 
In this case we have the additional relation, between t,(P) and the mean curvature 
H(P) of the surface of constant 7, at a given point P of the porous medium (Ericksen, 1963): 
V et,(P)=-2H(P) [11] 
From [10] and [11] it follows that: 


sg 


3 EO) ds 
V:(Q=v, (Pye * [12] 
For an isotropic medium, from equation [2] and [6] it follows that 

~ Ou. 

= [13] 
Then, from [12] and [13] we have: 
8Q 
2 ES) ds 

p,(Q) = @(Q) K(P) e* [14] 


p,(P) @(P) K(Q) 
To apply this last equation to study the erosion process in the filter, we need 
information about the filtration lines and the distribution of mean curvatures of the surfaces of 
equal z, along the filtration lines. 


3.3 Hydrostatic pressure and osmotic flow: thermodynamic arguments with a kinetic 
foundation. 


The chemical potential of water varies if the hydrostatic pressure or the water activity 
(proportional to water concentration in a solution) change (at constant temperature in our 
case). (The foundations of thermodynamic approach to nonequilibrium processes are given, 
for example, in Katchalsky and Curran, 1967). 

If we have pure water in the left chamber and water with disolved sucrose in the right 
chamber, but with the same hydrostatic pressure at the filter level, the chemical potential of 
water in the left side of the filter is higher than this potential in the right side, and it must 
decrease continuously, in the direction of the osmotic flow, through the interconnected pore 
space inside the filter. If there is no sucrose inside the filter, taking into account that water is 
scarcely compressible (so that water activity doesn’t change significantly in the pore space), 
the decrement in the chemical potential of water must be produced mainly by a fall in 
hydrostatic pressure inside the pore space. But as hydrostatic pressure is the same at the filter 
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level in both chambers, at the exit of the filter it must appear a sudden increment of 
hydrostatic pressure, equal and opposite to the fall in pressure produced during the filtration 
process. This counter-pressure opposes to the effects of the filtration pressure in the material 
located at the interface. Now, it is necessary to refine the argument in two directions. 

First, we don’t have pure water inside the filter, because a solution of 1mM KCl was 
used for the electroosmotic experiments, with the same concentration of KCl in both 
chambers and in equilibrium with the filter. As the surface of the grains was negatively 
charged, the entrance of K* was enhanced and that of Cl’ was hindered by the material of the 
matrix. This could be confirmed by the direction of the volume flow associated with an 
applied electric current through the filter. But the magnitud of the flow relative to the current 
suggested that the ions dragged not only their primary hydration cells, but their secondary and 
tertiary cells also, so that they were moving through a relatively wide pore space. In that case 
we can expect that the water potential is not significantly modified by the small scale 
fluctuations in ion concentrations inside the pore space, so that they don’t influence 
significantly the purely osmotic flow. 

Second, if we have a relatively wide pore space, even if an osmotic volume flow is 
pushing the substance out of the pores, sucrose can’t be totally excluded from the region of 
the filter near the interface. In this case the opposing hydrostatic pressure is established 
through a more or less narrow mixing width inside the filter. The numerical order of its length 
6 can be estimated by 

iis [15] 
(ee, ly 
where D, is an average of the dispersion coefficient of sucrose in the bulk of pore space, o, 
is the reflection coefficient of nonequilibrium thermodynamics and J, is the average volume 
flux density (equation [15] may be derived assuming an approximate equality between 
dispersion and advection of sucrose inside the pore space of the filter). 

To ground these thermodynamic arguments we need a relation between pressure, 
particle concentrations, temperature and forces between particles. Neglecting ionic effects and 
solid-liquid interaction at the boundary between the pore space and the matrix, and assuming 
that the forces can be expressed as a sum of forces between two particles, we have the 
following approximate equation (Green, 1969; Hobbie, 1974): 


pe ki(Cy+ C,) + [ few (0 ft 7 Cae 
0 


+ (2# fe. (r) f,,.@) 1 7 C.Cark (= fe. (r) f,, (x) r 7 Ge [16] 


Here k is Boltzmann’s constant; T is absolute temperature; C,, is the local number of water 
molecules per unit volume; C, is the local number of sucrose molecules per unit volume; 
Saw), Zys(r) and g..(r) are the normalized radial distribution functions of water molecules 


around a water molecule, of sucrose molecules around water molecules (or of water 
molecules around a sucrose molecule), and of sucrose molecules around a sucrose molecule, 
respectively; and f,.(r), f,,(@) and f,,(r) are the forces between two water molecules, 


between a water molecule and a sucrose molecule and between two sucrose molecules 
respectively. The radial distribution functions depend of T,C,,, and C,. 


If C, = 0, equation [16] would give us, in principle, an approximation to the field of 
hydrostatic pressure in the bulk of the filter, and when C, ¥ 0, it would give us the same 
field in the mixing width adjacent the sucrose chamber. 
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For pure water at 300° K (which is approximately the temperature of these 
experiments) and when the hydrostatic pressure p is of the order of 10° Pa (that is, 1 


atmosphere), the term kTC,, is positive and of the order of 10° Pa. As consequence of this, 
the other term must be negative, and also of the order of 10° Pa. Furthermore, the integral 


+00 
3 [eee (3 Cay Tt. r dr 
0 


is very sensitive to small changes in the radial distribution function; for example, those 
changes produced by small variations in the concentration C,, of water molecules. 


4. FLOW, FRACTURE, EROSION AND OSMOSIS IN A GRANULAR POROUS 
BODY. 


4.1 Geometrical analysis of filtration pressure. 


If we assume that the porous medium is homogeneous and isotropic at Darcy’s scale, 
equation [3] reduces to Laplace’s equation. For a cylindrical filter with boundary conditions 
[4] and [5], the filtration potential z,is a linear function of the coordinate along the axis of 
Ap 
L? 
where Lis the length of the filter and Ap is the fall in average hydrostatic pressure associated 


the cylinder and the filtration pressure is everywhere constant and of magnitude p, = 


with the given volume flow. 
In the osmotic case, the average force per unit volume, directed against the flow and 


due to the rise in pressure through the sucrose mixing width, is givenby py = @¢ = where 


6 can be estimated by equation [15]. But the circular end faces of the filter are not perfectly 
flat, even at Darcy’s scale, because when the filter is cutted from the original capillary tube 
(full of matrix), an irregular surface of pull-out fracture is formed in the matrix at these end 
faces. At the same time, the local porosity is increased there. As a consequence, some cavities 
and protuberances appear at these interfaces. When water is filtering, a surface of constant 7, 
follows the boundary of those cavities and protrusions, corresponding to a free water 
condition at the ends of the filter. 

Now, let us follow a filtration line in the flow direction, when the driving force is a 
difference in hydrostatic pressure between the chambers. If the line ends in a flat portion of 
the matrix-free water interface, there is not a significant variation in p, along the line, when 
we pass from a point P, in the bulk of the filter to a point Q, in the interface (fig. 2), because 

Soi 
the integral [H(s)ds that appears in formula [14] will be nearly zero. If the filtration line 
Sp 
ends in the bottom of a cavity, then when we pass from a point P, in the bulk of the porous 
medium to a point Q, in the interface, p, can increase significantly, because in a 
neighborhood of the interface the line will cross the surfaces of constant 7,in points where 
the mean curvature is positive and non-negligible (fig. 2). 
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Fig.2 Representation of iso-filtration potential surfaces and the corresponding orthogonal 
filtration lines 


When the depth of the cavity increases and the mean curvature at Q, increases also, the 

integral  H(s) ds will increase too. The extrapolation of the results of the analysis of several 
Sp2 

idealized models with plane symmetry (for example, an infinite groove with half-elliptical 

cross section in the surface of a half-space filled with a porous medium, homogeneous and 

isotropic; see Harr, 1991) suggest that p;(Q,) must be an increasing function of the product 


of a measure of the depth of the cavity or tunnel by the mean curvature at the bottom. When 
the cavity is a half-sphere, this product is equal to one. From the solution of Laplace’s 
equation for an idealized model (Wallace, 1972) we find that p,(Q,) is three times 


p,;(P, ) for a spherical cavity. But for lengthen cavities p,(Q,) may be one or two orders of 
magnitude greater than p,(P,). However, if the filtration line goes from a point P, in the 
bulk of the medium, to a point Q, in the upper part of the cavity or in the top of a protrusion, 
where the mean curvature of the surfaces of constant 7, is negative, p»(Q;) will be less than 


Py (P3) (fig. 2). 
4.2 Grain pull-out by filtration pressure. 


Let us consider a grain located at the interface. Let V, be its volume and A, its area. 


Then, the force due to filtration that tends to separate it from the rest of the matrix can be 
estimated as @V, ps where pr is the local value of the filtration pressure (at Darcy’s scale) 


and the correction factor @ relates this average with the microscopic field of forces there (at 
pore’s scale). The resisting force due to the cohesive interactions between the given grain and 
its surroundings can be estimated by SA, 7, where # is a correction factor which gives the 


fraction of the total area effectively engaged in the interaction with the other grains, and 7 is 
a mean shearing stress. Then, if Q is the point where the given grain is located, and Pis a 


point in the bulk of the filter, that belongs to the same filtration line asQ, from 


AD 
pr(P) = a) F, 
from equation [14] and from the equilibrium condition 
av, P-(Q) = BA, 7(Q) 


we obtain the following approximate expression for the shearing stress: 
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sQ 
(QF ¥t 1D or ap U7) 
BA, L 
We further assume that when 7 reaches a critical value 7,, the grain is separated and 
thrown into free water. When Ap increases from zero, a given critical value of 7 will be 
attained first at points where the product of the first three factors in equation [17] is a 


maximum: usually this will occur in points located at the bottom of the deepest initial cavities 
of the matrix-free water interface. 


4.3 Tunneling and its osmotic blockade. 


As a consequence of all this, when the difference in hydrostatic pressure between the 
water chambers reaches a critical value, a shear failure will be produced at the surface 
between grains, and one or several grains will be separated from the interface. One or several 
selfreinforcing filtration-erosion processes will be thus produced. Usually one of them will 
dominate the others, attracting progressively to its bassin most of the volume flow through the 
filter. The corresponding tunnel will grow faster and faster until it reaches the free water of 
the left chamber and the hydraulic short-circuit is completed. 

However, when the driving force of the flow is a difference in osmotic pressure, the 
filtration pressure at the bottom of the cavities will be counteracted by another volume force, 
Pop, aS was explained in 3.3. This volume force may be several orders of magnitude greater 


than the averagep,, because 6 << L, so that it can usually block the beginning of the 
erosion process in all the cavities. 


5.1 FLOW AND FRACTURE IN A POROUS BODY WITH A CONTINUOUS 
MATRIX. 


From equation [9] it follows that the tensile components of the matrix stress tensor 
will be larger in locations where the product @p is higher, if we assume static conditions 
(Fliigge, 1972). If we follow a given filtration line, the average pore pressure p will be a 
decreasing function of the arc length, so that if @ doesn’t increase too much and the other 
conditions remain constant, the tensile components of the matrix stress tensor will decrease 
along the filtration line. Now, consider a wide master crack in a saturated porous beam with a 
continuous matrix, such that the crack reaches the lower boundary of the body. If the beam is 
in the air and is flexed as suggested in fig. 3, water will tend to flow downwards. 


MASTER CRACK 


Fig. 3 Filtration lines and surfaces of constant filtration potential for a flexed beam with a 
master crack in its tensile field. 


The water inside the wide crack can be considered as free water. Then a surface of constant 
filtration potential will follow closely the lateral faces of the crack, and the other surfaces of 
constant 7, will be disposed as is shown qualitatively in the figure. The tip of the crack is ina 


tensile nominal stress field due to the external loads. However the magnitude of the tensile 
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components will be disminished by the effect of the filtration forces (that tend to close the 
crack), due to the decrease of p along the filtration lines. So, the flexural strength of the 
porous beam will be higher when it is wet, relative to its value when the body is dry. This 
difference should increase when the average porosity increases. 

To study the compressive strength, let us consider a vertical cylinder of saturated 
porous body, with a population of spherical cavities (Fig.4). 


j ® | eee 
1% 7 he oo ae 


A. gis | Stress 
microcrack 


fF nominal - 
‘tensile 


Fig. 4 Compressed cilinder and stress distribution around a cavity with microcracks. 


The radii of these cavities are at least of the same order that the Darcy’s scale for the porous 
medium, so that they do not belong to the interconnected pore space used to construct the 
averaged fields. Many little cracks begin at the surfaces of the cavities, and go into the matrix, 
in all directions. If we apply compresive external loads to the flat faces, a radial filtration 
pressure field will be produced, directed toward the curved face which is both free from stress 
and a free-water boundary of the porous body. Then, the nominal stress field is the 
superposition of a vertical compressive stress (due to external load) and a horizontal radial 
tensile stress field (due to filtration). The compressive stress by itself produces a tensile stress 
in the matrix surrounding the upper and lower parts of the spherical cavities. If this stress is 
big enough, a crack located there will grow and a fracture process will begin. The radial 
tensile field will produce an increment of the tensile stress already induced in the matrix by 
the compressive external load. As consequence, due to flitration, the compressive strength of 
the wet porous body will be less than its compressive strength when it is dry. 


6. CONCLUSIONS. 


(a) The pull-out of the grains due to the filtration force at the interface, the subsequent self- 
reinforcing growth of an erosion tunnel in the bulk of the filter, until a true hydraulic short- 
circuit is produced when the volume flow is of a purely hydraulic origin, as well as the 
blockade of this process from its very beginnings, due to a counter-force produced in the same 
interface when the flow is due to an osmotic thermodynamic force, seem to be a complete and 
coherent explanation of all the observed experimental facts. However, new experiments 
and/or digital simulations could be suitably designed and done to isolate and separately 
confirm the main parts of the argument. 

(b) It is known that the flexural strength of a wet cement paste is usually higher than the 
flexural strength of this same cement paste after it is dried. The opposite happens with the 
compressive strength (Neville, 1970). Taking into account the differences between the failure 
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mechanisms of cement paste under tension and its failure mechanisms under compression 
(Wittmann, 1993), the analysis presented in section 5 could be used as a starting point to 
develop a quantitative theoretical approach to the relation between water flow and fracture in 
saturated porous concrete. 

(c) In principle, from equation [16] we can give a molecular foundation for the assumed fall in 
hydrostatic pressure inside the interconnected pore space of the filter, and also to the rise in 
pressure through the mixing width when there is sucrose in the right chamber. However, a full 
development of this molecular approach requires knowledge of the radial distribution 
functions and of the fields of concentration of water and of sucrose in the pores. By now, this 
approach does not seem to be feasible. A less detailed approach would be to construct an 
averaged version of equation [16], as well as averaged versions of the transport equations for 
water and sucrose, always at Darcy’s scale. 
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ABSTRACT: The fracture process of a continuous matrix in a porous medium under the combined 
effect of filtration and external mechanical loads is considered. Taking into account the differences 
between the failure mechanisms of cement paste under tension and its failure mechanisms under 
compression, an analytical approach to the relation between water flow and fracture in saturated 
porous Portland cement pastes is developed. The well known differences in behaviour between the 
flexural and compressive strengths of wet and dry Portland cement pastes is explained. The exten- 
sion of the obtained results to the flexural and compressive strength of normal concrete is briefly 
discussed, including suggestions for further experimental and digital simulation work 


1 INTRODUCTION 


Hardened Portland cement paste is a porous and viscoelastic material whose anelastic 
properties depend strongly of the relative humidity and water redistribution processes in 
the gel pores. The smaller gel pores (mainly of sub-micrometer dimensions) and the 
technological pores (from small capillary pores with dimensions of the order of 10 4am 
to big crack-like pores with lengths of the order of 1mmor more) form a system of 
communicating cavities or channels, bifurcating and recombining at random. Under ap- 
propriate conditions, water flows through these tiny and irregular channels, dominated 
by viscous friction. Flowing water exerts both normal and shear stresses on the surface 
of the solid matrix, like an external load. After averaging to smooth the microscopic 
spatial fluctuations of the fields, this mechanical effect can be described by the so called 
“filtration pressure”, a force per unit volume of the porous medium, acting on the solid 
matrix and produced by the flow. If there are closed cavities (like trapped air bubbles) 
that remain outside the flow, they must be considered together with the solid matrix, 
although they modify the local stress and strain distributions. 

Thus, the strain and stresses produced by the filtration process must be combined with 
the strain and stresses produced by the external loads to describe the mechanical failure 
of a wet porous medium, like saturated cement paste, mortar or concrete. 


It is well known that the flexural strength of a wet cement paste is usually higher than 
the flexural strength of this same cement paste after it is dried. The opposite happens 
with the compressive strength (Neville, 1995). 


Considering the differences between the failure mechanisms of cement paste under ten- 
sion and its failure mechanisms under compression, the analysis of the fracture process 
in a porous medium is used here as a starting point to develop an analytical approach to 
the relation between water flow and fracture in saturated porous cement pastes. 


This approach is then applied to study the relation between the field of filtration pres- 
sures and: 


(a) The compressive strength of cylinders of cement paste, mortar or concrete (used as 
standard test bodies). 


(b) The flexural strength of beams of these same materials (that are also used as stand- 
ard test bodies in four-point bending). 


2 MODELLING TOOLS 


To model the aforesaid filtration and fracture process in the cement paste we will use 
three kinds of tools: 


(1) An average form of the equations that govern the transport processes in the pore 
space and the stresses in the solid matrix (Hall and Hoff, 2002; Bear, 1972; 
Scheidegger, 1963). 


(2) Some results for solenoidal and almost gradient fields, taken from classical field 
theory and differential geometry, and applied to the average field variables (Ericksen, 
1960; Suarez-Antola, 1983; Suarez-Antola, 1997). 


(3) A certain population of defects included in the matrix and a quantitative expression 
of the commonly admitted mechanisms of fracture of cement pastes under flexural and 
compressive loads (Chtchourov, 1980; Wittmann, 1983). 


2.1 Average field equations 


The fields of velocity and stress in the fluid and the stress field in the solid matrix fluc- 
tuate randomly along distances of the order of the local pore’s diameter d (for a defini- 
tion of the local pore diameter see Scheidegger, 1963). 

However, it is possible to smooth out these fluctuations by averaging the field variables 
over regions of characteristic length7 , much greater than d and called Darcy’s scale. 
Selecting /, at its turn, much less than a characteristic length Lof the whole porous 
body, we obtain averages that can be considered as true local values. 

A possible choice for / is the geometric mean / = Vd.L.. 

In our case we can estimate d=10° m or less, andL =10°! m or less, so that /=JdL 
between 10“*and 10° m may be a suitable length scale for averaging purposes. 

For each point in the porous medium we consider an averaging region that has this point 
as its centre, and then we assign to the point the averaged values of all the field varia- 
bles, with independence of the true location of the point (in interconnected porous 
space, in solid matrix, including isolated voids, or in the boundary between them). 


We thus obtain a vector field of filtration velocities V,;, a scalar field of averaged pore 
hydrostatic pressure p, a scalar field of effective porosity ¢ and a tensor field of aver- 
aged matrix stressesc . 

The emergent equations at Darcy scale are: 


a __v. 
ea (1) 
eee (2) 
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where Ve() represents the divergence, V() the gradient, ethe scalar product, 
z>=pt+pgz is the filtration potential ( is water density, g is gravitational accelera- 
tion, z is the vertical position from a reference horizontal plane, z increasing upwards), 
Lt is water’s dynamic viscosity and K is a scalar field that express the local hydraulic 
conductivity of the porous medium . 

Equation (2) gives Darcy’s law for an isotropic medium. 

Equation (1) connects the variation of the fraction of pore’s volume with water flow 
(Jaeger, 1969). 

If ois the tensor of averaged stresses in the matrix, the force acting on an element of 
surface of area AA and unit normal vectorn, at the scale of Darcy is given by the for- 
mula: 


AF=AA(6-pdlDen (3) 
Then the complete average stress tensor in the solid matrix verifies the emergent static 
equation Ve(o-gpl)=X (4) 


Here I is the unit tensor and X is the average external force per unit volume. The av- 
erage force on the matrix, per unit of volume of a saturated porous medium, due to wa- 
ter flow, is the filtration pressure: 

Pp =-9V 2, (5) 
So pis everywhere orthogonal to a surface of constant z , . The filtration pressure is re- 
lated with the average matrix stresses due to the filtration process o ; by: 

p co 7 Ve on (6) 


2.2 Solenoidal and almost gradient fields 


For an isotropic medium, from equation (2) and (5) it follows that 


Py =fh4, (7) 


If the phenomena of interest occur in time scales that are much less than the time scale 
of variation of pore’s volume, equation (1) reduces to this one: 

Vev, =0 (8) 
Let us consider a vector line of the fieldV,. If t; (s) is the unit tangent vector at the 
point of intrinsic coordinate s (s being the arc length), given two points Q and P that 


belong to the line, we may apply Bjorgum’s characterization of solenoidal (divergence- 
free) fields (Ericksen, 1960; Suarez-Antola, 1983, Suarez-Antola, 1997) 


sQ 
- vet (s) ds 


V/(QQ)=v,(P)e * (9) 


(Here v, =|/v,|| is the magnitude of the vector field). 
Now, according to equation (2), V,is an almost gradient field. 


Then each tangent unit vector t, at a point of a filtration line is at the same time a unit 
normal vector to the surface of constant filtration potential that pass through the same 
point. r 

In this case we have the additional relation, between t,(P) and the mean curvature 
H(P) of the surface of constant 7, at a given point P of the porous medium: 


Vet,(P)=-2H(?) (10) 
From (9) and (10) it follows that: 
2 fi) ds 
Vs(Q)=v,(P)e * (11) 


sQ 


PQ _ sk) 7)" 
p-(P) ¢(P) K(Q) 
To apply this last equation to study the fracture process in the cement paste, we need 


information about the filtration lines and the distribution of mean curvatures of the sur- 
faces of constant z, distributed along and orthogonal to the filtration lines. 


Then, from (7) and (11) we have: 


(12) 


From equations (2) and (8) it follows that the filtration potential verifies the elliptical 
equation: Ve (KVz, )= 0 (13) 


At an impermeable boundary (as the interface between the steel plates of the compres- 
sion machine and the tested standard cylinder or cube) Vv, en =0, ni being the unit nor- 
mal vector to the interface at the considered point, so 

ne Vz, =0 (14a) 


At free water interface (as the free interfaces of the cylinders, cubes or beams used in 
standard strength tests: 
7, =cons tant (14b) 


Equation (13) jointly with the boundary conditions (14a) and (14b) for z,, allows us to 
determine the surfaces that correspond to constant values of 7, . 


2.3. Population of matrix defects and fracture under flexural and compressive loads 


The characteristic length / of Darcy’s scale may be used to classify the population of de- 
fects in two broad classes: the microscopic defects whose sizes are much smaller 
than/ and the macroscopic ones, much larger than/ . 

The usually small fraction of defects with the same order of magnitude than / may be 
called mesoscopic. 

At a microscopic level the irregular fracture surface of the cement paste seems to ad- 
vance through a pull-out process that separates the lenticular crystals of hydrated calci- 
um silicate from each other, breaking the van der Waals bonds between these fairly 
elastic elements. 

But from a macroscopic point of view the failure under simple compression begins with 
the growing of many micro-cracks from small pores, parallel to the external load, and 
continues with the fusion of the pores in master cracks (Zaitsev, 1980). 

Thus, the compressive strength depends strongly of the cement paste porosity. 

Also from a macroscopic standpoint, the failure under tension often begins with the 
growth of a master crack already located in a highly stressed volume (HSV, having 


tensile stresses between 95 and 100% of the maximum) and often continues with rela- 
tively few accompanying small cracks. 

So, the tensile strength depends weakly of the cement paste porosity, and strongly of the 
lengths of the big crack-like pores in the HSV (Kendall et al., 1983). 

Size and shape of the body, the load distribution applied to it and other systemic proper- 
ties (Bazant, 2004), have a strong influence on the fracture process and the measured 
strength values. 

So, two specific cases of practical interest will be considered first: the test of the com- 
pressive strength of standard cylinders and the test of flexural strength of standard 
beams of square cross section. 


3 SOME SIMPLE MODELS OF FLOW AND FRACTURE IN CYLINDERS 
AND BEAMS OF CEMENT PASTE 


3.1 Filtration effects in uniaxial compression of a porous cylinder 


To study the compressive strength, let us consider a vertical cylinder of saturated porous 
body, with a population of spherical cavities (Figure 1). 


+} 


Figure1. Sketch of a compressed porous cylinder with a population of spherical pores 


The radii of these cavities are at least of the same order that the Darcy’s scale for the po- 
rous medium, so that they don’t belong to the interconnected pore space used to con- 
struct the averaged fields. 

Many little cracks begin at the surfaces of the cavities, and go into the matrix, in all di- 
rections. If we apply compressive external loads to the flat faces, a radial filtration 
pressure field will be produced, directed toward the curved face which is both free from 
stress and a free-water boundary of the porous body. 

The vertical component of filtration field, due to gravity, may be neglected in the mod- 
elling of fracture in compressive strength tests. 

Then, the nominal stress field is the superposition of a vertical compressive stress (due 
to external load) and a horizontal plane tensile stress field (due to filtration) with both 
radialo ,,, and tangentialo,,, principal components. 


The compressive stress o, by itself produces a tensile stress in the matrix surrounding 
the upper and lower parts of the spherical cavities. 

If this stress is big enough, a crack located there will grow and a fracture process will 
begin. 

The radial and tangential tensile fields will produce an increment of the tensile stress al- 
ready induced in the matrix by the compressive external load. 


As a consequence, due to filtration, the compressive strength of the wet porous body 
will be less than its compressive strength when it is dry. 

The symmetry of the problem simplifies the derivation of formulae for the filtration 
stresses o ,,, anda, , once the field of filtration pressures p, (t,r) is known. 


Let 7, be the radius and / the height of the standard cylinder. 
Applying a volume balance to a coaxial cylinder of radiusr less thany,, and heights, 


(2achv, =—-ar’.h. ) we obtain the field of radial filtration velocities 
t 
r 0o 
_=- =, 15 
i ear (15) 
From (15) and (7), it follows that 
Py (t, r) = f(t)r (16a) 
Up OP 
t)=-—— .— 16b 
f\ ) 2.K Ot oe 


As the time scale of the filtration process is at least an order of magnitude greater than 
the duration of the strength test, during the period of load application f (t) xf (0). 


Om, _ HM Op 
Or 2.K Ot 


But from (5) it follows that: 


Integrating between 0 andy, and taking into account that 7, (t, Yo ) © Ty (0, 1 ) =, 18 


the atmospheric pressure (neglecting gravity) and is thus negligible relative to 


Ty (1,0) = I y (0,0) = lo,|, we obtain: 
ry, Od 
|--Fo. (17) 
4K ot 
Then, from (16b) and (17), 
2 
f= 10)==4. (18) 
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Now, the problem of finding the filtration stresses when the filtration pressure verifies 
(16) and the solid matrix is homogeneous, isotropic and linear elastic is mathematically 
the same as the well known problem of finding the stresses in a linear elastic rotating 
disc (Nadeau, 1964). 


Adapting the solution of this classical problem to the filtration case we obtain (being 


v Poisson’s ratio): 
2 
wat |e P-(2] | (19a) 
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Figure 2 Sketch of a spherical cavity with a vertical crack of length c 


A vertical crack extends from the upper surface of the cavity, with its plane orthogonal 
to the radial direction (tangential crack). 

Let us suppose that the onset of crack propagation occurs when the nominal tensile 
stress field at the crack location reaches a threshold value c,, . 


If there is no filtration stress, the only nominal tensile stress is due to the compressive 


Cc 


2 


load Pp, and its value is |o,|=—S— 
ro 


as 


At the critical state, the threshold load is z.r,’.o,, = oh 


If there are filtration stresses, o, , has no effect in this case, but o,,. produces a nominal 


tensile stress 3.0,, at the location of the crack, so that now the total nominal tensile 


stress is|o, 


DG fos 


Taking into account equation (19a), the critical stress would be reached when 


143040) 4 1 [=] =, 


Yo 


lc. 


The new threshold load 


P| will be related with the old one by the relation: 


(20a) 


If the plane of the vertical crack is radial, o ,, has no effect and o , ,is active now. In 
this case we obtain from (19b): 


(20b) 


3.2 Filtration effects in bending of a porous beam in a flexural strength test 


Now, consider a wide master crack in a saturated porous beam with a continuous ma- 
trix, such that the crack reaches the lower boundary of the body. If the beam is in the air 
and is flexed as suggested in fig. 3, water will tend to flow downwards. 


Figure 3. Sketch of a flexed porous beam with some filtration lines (represented by continuous 
curves) and several surfaces of constant filtration potential (represented by dotted curves). 


When water is filtering, a congruence of regular surfaces of constant 7, fills the porous 


body. The flow velocity and the vector of filtration pressure are at every point orthogo- 
nal to the surface of constant filtration potential that passes through this point. At each 
boundary with free water conditions, a surface of constant 7, closely follows the 


boundary. The water inside the wide crack can be considered as free water. Then a sur- 
face of constant filtration potential will follow closely the lateral faces of the crack, and 
the other surfaces of constant z, will be disposed as is shown qualitatively in the figure. 


The tip of the crack is in a tensile nominal stress field due to the external loads. 


However the filtration lines will arrive orthogonally to the crack surface, so that the 
magnitude of the tensile components will be diminished by the effect of the filtration 
forces (that tend to close the crack). 

So, while the presence of the crack lowers the local elastic energy, the filtration pressure 
around the crack raises this local elastic energy and modifies the global energy balance 
that appears in the energy theory of fracture. 

Indeed, from equation (3) it follows that the tensile components of the matrix stress ten- 
sor will be larger in locations where the product ¢p is higher, if we assume static con- 


ditions. 


If we follow a given filtration line, the average pore pressure p will be a decreasing 
function of the arc length, so that if ¢ doesn’t increase too much and the other condi- 
tions remain constant, the tensile components of the matrix stress tensor will decrease 
along the filtration line. 

Beginning with equation (4), it is possible to derive an approximate formula for the 
elastic energy release in the solid matrix, due to the presence of the crack, including the 
effect of the filtration process. 

If the master crack crosses completely the beam of widtha, and if cis its length, then 
the elastic energy released is 


ae 
7 nerafs = ‘Fe. (21) 


Here o,, is a horizontal nominal tensile stress in the matrix, pis a nominal pressure in 


the pore space, and E is Young’s modulus. 

If G, is the average energy absorption in the fracture process zone (FPZ) per unit area 
of crack advance, the increase in the beam’s energy due to the presence of the crack 
is2.G,.a.c (22). 

The sum of these two energies will have a maximum when the stress in the matrix veri- 
fies the threshold condition: 


2.G,.E 
ox -2.¢.p.c0, -——— =0 (22) 
H.C 


Botho and pare proportional to the bending moment M produced by the loads dur- 
ing the strength test: 


o, =k,(z)M (23a) 
p=k,(t,z)M (23b) 
As the master crack length may be several mm, while the height a of the beam is not 


less than50mm ,, we may suppose that k, ~ a (J is the moment of area of the cross 


section). 


The dependence of ‘, on the system’s parameters remains to be studied. Nevertheless, 


for the same reasons that in the compressive strength test, we can suppose that 
k ,(t,z)=k,(0,z) 


From (22) and (23) we obtain: 


MM, »——“— (24) 


[2.G,.E , 
By definition: M,, = a ando,,, * ,,——— 1s a well-known stress threshold of frac- 
TC 


ture mechanics. As a consequence, the flexural strength of the porous beam will be 
higher when it is wet and saturated, relative to its value when the body is dry. This dif- 
ference should increase when the average porosity (of the pore space in which water 
flows) increases. If as a coarse approximation we putk, ~k,,, from (24) it follows that 


the quotient of threshold bending moments is only a function of the porosity of the pore 
space engaged in filtration. 


To quantify the concentration of filtration pressures due to the crack, let us follow a fil- 

tration line in the flow direction. If the line ends in a flat portion of the lower boundary 

of the beam, there is not a significant variation in p, along the line, when we pass from 

a point P, in the bulk of the medium to a point Q, in the interface (fig. 3), because the 
Sai 

integral {H(s) ds that appears in formula (12) will be nearly zero. 


Sp] 


If the filtration line ends in the bottom of a cavity, then when we pass from a point P, in 

the bulk of the porous medium to a point Q, in the interface, p, can increase signifi- 

cantly, because in a neighbourhood of the interface the line will cross the surfaces of 

constant 7 ,in points where the mean curvature is positive and non-negligible (fig. 3). 

When the depth of the cavity increases and the mean curvature at Q, increases also, the 
Soo 

integral [H(s) ds will increase too. 


Sp2 


The extrapolation of the results of the analysis of several idealized models with plane 
symmetry (for example, an infinite groove with half-elliptical cross section in the sur- 
face of a half-space filled with a porous medium, homogeneous and isotropic; see Harr, 
1991) suggest that p-(Q,) must be an increasing function of the product of a measure 
of the depth of the groove (or notch) by the mean curvature at the bottom. For lengthen 
cavities like cracks, p;(Q,) may be one or two orders of magnitude greater 
than p;(P,). 


4 MORTARS AND PLAIN CONCRETES 


Up to now, we studied the onset of crack propagation under the combined effects of ex- 
ternal loads and filtration pressure, in two idealized situations: a spherical cavity with a 
vertical crack in a compressed cylinder of hardened cement paste, and a vertical master 
crack in a beam of cement paste in four points bending. 

In order to extend the analysis of the effects of the filtration pressure to the cases of 
compressive and tensile failures of normal mortars and concretes, several things must be 
considered. 


Amongst them: 


(1) The transition ring, that appears at the interface between the paste and the aggre- 
gates (Maso, 1980). 


(2) The difference in stiffness between the cement paste and the fine and coarse ag- 
gregates. 


(3) The gradual strain softening and the crack band in the FPZ (Bazant and Oh, 
1983). 


(4) The stages of crack growth prior to fracture: initiation, slow stable crack growth, 
crack arrest, a true threshold condition, and unstable crack propagation up to definitive 
rupture. 


Let us consider the typical case of a mortar or concrete with relatively non-porous ag- 
gregates, stiffer and stronger than the cement paste that surrounds them. 

At the time of mixing, a film of water appears between the surface of the non-porous 
aggregate and the bulk of the water-cement mix. 

Due to this circumstance a transition ring of cement paste is formed, usually weaker and 
more porous than the bulk of the paste. 


If the aggregates are stiffer than the bulk of the cement paste, this usually produces 
stress concentrations in the aggregates and in the adjacent transition rings. In a nominal 
one dimensional stress field, like the one found in a standard compression test of a cyl- 
inder of mortar or plain concrete (neglecting the complications introduced by the restric- 
tive action of the plates of the testing machine), or in a small enough region of a beam 
in the standard four-point bending test, a non-uniform three dimensional state of stress 
is produced in the composite material due to the abovementioned difference between the 
elastic moduli of the paste and the aggregates. 


If the aggregates are stronger than the hardened cement paste, a growing crack will run 
around the aggregate travelling in the transition ring and eventually leaving it behind if 
the crack growth continues. Then, the crack path in normal concretes will be tortuous, 
and the aggregate particles will act as crack arresters. In high strength concretes things 
are usually different (Wittmann, 2002). 


In both, mortars and normal concretes, the fracture process begins with the growing of 
one or more cracks suitably located and oriented in relation with the stress field. Satel- 
lite micro- cracks emerge from the growing tortuous crack. A band-like fracture process 
zone (associated with the strain softening already mentioned) is produced. Because of 
the crack arresting processes and the complex pattern of cracks that is produced, hard- 
ened mortars and concretes are usually less brittle than hardened cement pastes. 
Concretes show two critical loads in standard compressive strength tests. 

The first critical load corresponds to crack growth initiation. In concrete, above the first 
critical load, Poisson’s modulus begins to increase monotonically towards 0.5. 

Between the first and the second critical loads there is stable crack growth. At the sec- 
ond critical load the cubic dilatation of the stressed cylinder reaches its maximum value. 
Above the second critical load there is unstable crack propagation to definitive rupture. 
However, the effect of the filtration field in the fracture processes of mortars and con- 
cretes, at least up to the first critical load, should be in its main traits fairly similar to the 
effect already discussed for cement pastes, not only for the compression test, but also 
(and mainly) for the bending test. 


Consider first a vertical cylinder of concrete, porous and wet saturated, under uniaxial 
compression. Due to the assumed differences in stiffness, the bulk of cement paste is 
less stressed than the transition rings, and these are less stressed than the cores of the 


aggregate particles. These variations are mild in relation with sand particles, but very 
significant in the case of coarse aggregates. 

The fields of filtration velocities and filtration pressures are locally modified by the 
non-uniformities in the stresses (that produce variations in the dimensions of the inter- 
connected pore space), by the stress-independent variations in porosity (higher in the 
transition rings than in the bulk of the matrix), and by the blocking effect of the imper- 
vious aggregate particles over the water flow in the pores. 

Nevertheless, on the average, they point outwards from the cylinder axis, as in the case 
of the cement paste, and tend to enhance the growth of cracks, mainly in vertical direc- 
tion. Thus, the first critical load in compression must be smaller for a wet cylinder in 
comparison with a dry one (all the other things remaining equal). 

We could expect that for a splitting mode of crack growth, the difference in the first 
critical load should increase with porosity in a way similar to the cement paste threshold 
load studied above. Between the first and the second critical load, the changing and ex- 
panding pattern of cracks and micro-cracks produces a progressive modification of the 
filtration field and the corresponding filtration pressures. 


The filtration field and filtration pressures in bended beams of mortars and concretes are 
locally modified by the heterogeneities as were described for a cylinder under compres- 
sion. 

However, the lower portion of the beam, where the fracture process begins with the 
growing of a master crack, is under tensile nominal stress. 

On the average this enhances the filtration process and increases the filtration pressures 
that tend to close the crack for the same reasons already considered for the cement paste 
beam. 

Nevertheless, the arresting effect of the aggregates and the transition rings increases the 
number and tortuosity of the pattern of satellite micro-cracks connected with the master 
crack. In spite of the blocking effect of the impervious aggregates, this could enhance 
the local filtration process even more. A 

s a consequence, for a wet saturated beam of mortar or concrete, the critical value of the 
bending moment that corresponds to the onset of a master crack growth should be great- 
er than the critical value for a dry beam, all the other things remaining equal. 


CONCLUSIONS 


Now we have a solution to the problem posed in the introduction: to explain the differ- 
ences in the behaviour of the flexural and the compressive strengths of a body of cement 
paste when it is tested wet and saturated or when it is tested dry. 

In the case of compressed cylinders, given a crack with a tangential plane and a thresh- 
old stress, from (20a) it follows that the most dangerous cracks from the standpoint of 
the onset of crack propagation are the ones located on the axis of the cylinder. 

The tangential cracks located near the surface of the cylinder are almost not affected by 
filtration stresses. Given a crack with a radial plane, from (20b) it follows that again the 
most dangerous cracks are located on the axis. 

However, now the cracks located near the surface of the cylinder are affected by filtra- 
tion stresses, and if they size is big enough, a fracture onset may happen in one of these 
radial cracks. In the case of bended beams, the model should be completed with a for- 
mula fork, . 

In principle this could be done from a water balance as was done for the cylinder, but 
the arguments are subtler and lengthier and will be developed elsewhere. 


Albeit the analytical formulae obtained in this work are crude approximations to very 
complex situations, they can be improved by means of more realistic analytical models 
of the onset of crack propagation. 

The energetic approach used here could be extended to cope with the essentially hetero- 
geneous nature of fracture processes in cement pastes, mortars and concretes. 

The obtained results may be compared with true experimental results, or with the results 
of digital simulations done with realistic nonlinear models and suitable computer codes. 


Even in the case of cylinders of cement paste, compressed without filtration pressure ef- 
fects, it is well known that the pattern of crack growth that leads to fracture may some- 
times correspond to what seems to be a shear mode. 

In fact, if the pores are randomly distributed, the pore cracks are randomly oriented and 
the pores are near enough, both digital simulations and experiments suggest that the 
probability of having a pattern of crack propagation that ultimately leads to rupture in an 
apparent shear mode instead of a vertical splitting mode, increases with the pore density 
(Zaitsev, 1980). 

From the qualitative discussion of the fracture processes of mortars and normal con- 
cretes under filtration pressure fields, given in part 4, it follows that the analytical ap- 
proach could be only a guide to further experimental and numerical studies in the sub- 
ject of the present paper. 

Only a thorough numerical approach can afford suitable tools and realistic results for an 
in-depth study of the fracture processes under the effects of filtration pressure fields in 
concrete, especially in problems related with structures. 
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